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1. Introduction 



The heavy quark is approximated by a static propagator, while for the light quark propagator 
we use estimated all-to-all propagators. To improve these estimates we use the so-called Domain 
Decomposition Improvement [1]. In these proceedings we describe only briefly the methods which 
we use to extract mass splittings, ratios of decay constant, and kinetic corrections for static-light 
hadrons. Also we show only a subset of our results. A complete discussion of the methods and our 
results is presented in Ref. [2]. 



2. Methodology 

2.1 Masses 

The masses of ground and excited states are extracted with the variational method [3]. For 
that purpose we construct the different interpolating fields for the states by using different "wave- 
functions" for the light-quark source and sink. The cross-correlation matrix C(t), which we can 
compute from them, is then inserted into the generalized eigenvalue problem: 



C{t) Y {a) =X ia) {t,t )C{t ) Y {a) ■ 
For sufficiently large t S> to, the eigenvalues are 

X^(t,t )=c^e-^ Mia) \l + 0(e-^ Aia) ) 



(2.1) 



(2.2) 



where is the energy difference to the closest state. However, due to the static quark propagator 
only differences to the so extracted lattice energies are physically meaningful. Thus, we only 
report mass splittings in the following. 

2.2 Couplings 

Using the eigenvectors obtained by solving Eq. (2.1), we can construct and fit the ratio 
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(2.3) 



which allows us to determine the couplings v. Ratios of different couplings to the same mass 
eigenstate [4] are even more straightforward: 
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Then, for example, the coupling of the local vector operator (0, = qjiQ, where q and Q denote the 
light and heavy quark, respectively) can be related to the pseudoscalar decay constants via 



PS 



(2.5) 
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In order to cancel renormalization constants and matching coefficients between HQET and QCD, 
we deal with the ratios of of decay constants. So, for example, the ratio f B > /fs s , may be extracted 
from the m q = m s point of 



f (2) (2) 
J PS _ V i 



M,, 



B[ ' +ff(ki 2)2 /m b ,ki l)2 /m b \. (2.6) 



where M w = 5400 MeV and the vj a) and E^ come from fits to Eq. (2.3) for a = 1 and 2. 

2.3 Kinetic corrections 

The G{\jmQ) kinetic corrections to the static approximation can be incorporated into the 
simulations in form of lattice three-point functions 

T(t,t')ij = (0\qOiQ(t) QD 2 Q{t') QO]q(0)\0) , (2.7) 

with the current insertion QD 2 Q(t'), where D 2 is the lattice-discretized covariant Laplacian. To 
obtain the corrections not only for the lowest lying state but also for the excitations, we consider 
two separate variational problems: 

C(t - t')ij = ( 1 qOiQ{t) QO)q{t') 1 ) , 

C{t')ij = {Q\qOiQ{t') QO)q{0)\0). (2.8) 

Solving the two corresponding generalized eigenvalue equations, 

C(t - 1') =X(t-t'4- 1')^ C(t' - 1') y {a) , 

C(t')^=X(t',to) {p) C(t )f P) . (2.9) 

then gives sets of eigenvectors, which can be used to project the states of interest in Eq. (2.7). To 
cancel exponentials and some overlap factors one may form ratios: 

where the e( a 'W represent the matrix elements relevant for the kinetic corrections. 

3. Simulation details 

Our calculations are performed using CI fermions [5] for the light valence quarks. For the 
quenched simulations we work with three lattice sizes with approximately the same spatial volume 
of 2.4 fm and lattice spacings reaching from 0.20-0. 12 fm. A necessary step for the CI operator is to 
smear the configurations. For the quenched lattices one step of HYP-smearing is applied. Another 
part of our simulations are done on two sets of lattices with Nf = 2 dynamical CTfermions [6]. 
Here one step of stout smearing is used. In all cases we use the one-loop improved Luscher-Weisz 
gauge action [7]. 

In order to increase the number of basis operator for the variational method, we create quark 
sources and sinks with a number "shapes" by using different covariant spatial smearings. 
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Figure 1: Meson and baryon masses, relative to the meson ground state M\$ = M(B q '), as a function of M\ 
(oc m q ) on the 16 3 x 32, j3 = 7.90 quenched lattice (left plot) and the 16 3 x 32, j3 = 4.65 dynamical lattice 
(right plot). All masses result from fits to the eigenvalues of the complete 4x4 basis. Circles represent the 
S states; squares, the P states; diamonds, the P + states; bursts, the ID±; pluses, the baryons Ag and Eg ; 
and crosses, the experimental results [8] (using ro = 0.49 fm). 
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Figure 2: Quenched continuum extrapolations of excited B s meson masses, relative to B s (left plot) and 
quenched continuum extrapolations of the Ei, —A/, and Q.Y —Kb baryon mass differences (right plot). 



4. Results for mass splittings 

In the following we present a selection of our results for the static-light hadrons mass splittings. 
Figure 1 shows exemplarily the results obtained on the 16 3 x 32, jS = 7.90 quenched lattice (left 
plot) and the 16 3 x 32, j8 = 4.65 dynamical lattice (right plot). The vertical dashed lines indicate 
the physical pion mass and the pion mass which corresponds to the strange quark mass on these 
lattices. The latter is set via the splitting M\$ s — M\s ud = 76.9 MeV, which is the 1 /M H (*) — > linear 
extrapolation of the experimental values M w — M g (») = 86.8 MeV and M m — M D («) = 103.5 MeV. 
On both lattices we can extract a number of excited states, including a 3S state in the quenched case. 
The results on the dynamical lattices are consistent with those on the quenched configurations but 
have much larger error. This is most likely due to the additional fluctuations which are added by 
the sea quarks. 

On the left hand side of Figure 2, we show continuum extrapolations for our quenched results 
for mass splitting of the B s states. While the results on the two finer lattice agree with each other we 
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Table 1: Summary of mass differences (in MeV). When possible, we compare our results to experimental 
numbers [8]. Our values are given using ro = 0.49( 1) fm. The first error of our results is the statistical error 
while the second one is a systematic error coming from the uncertainty of ro. 



P (jB,/f B ) static static {fir, I fa) static 



7.57 


1.042(7) 


1.309(75) 


0.976(142) 


7.90 


1.058(21) 


1.237(15) 


0.996(59) 


8.15 


1.075(26) 


1.259(25) 


0.972(61) 


oo 


1.087(31) 


1.240(58) 


0.972(123) 


4.65 


1.108(29) 


1.356(142) 


1.089(259) 


5.20 


1.089(41) 


1.453(168) 


1.026(128) 



Table 2: Static-light decay constant ratios and the quenched continuum values. Values are given using 
ro = 0.49(1) fm to set the physical m s point. 

find a significant difference for the coarsest lattice, especially for the 2S state. This might be a hint 
for possible discretization errors on that lattice. Nevertheless, for the continuum extrapolation we 
include also the results on the coarsest lattice. The right plot in Figure 2 shows the corresponding 
extrapolation for mass differences of some of the baryons alone. 

Table 1 summarizes the results of our simulations in physical units and compares (where pos- 
sible) with experimental values. 

5. Results for decay constants 

In this section we present first results for the ratios of meson decay constants using the method 
described in Section 2.2. Figure 3 shows our results for the ratios of meson decay constants 
(/b s //s) static (left plot) and (Jb'J Ib s ) static (right plot) as a function of lattice spacing. For the 
continuum extrapolation of the quenched results we try different fit functions. Apart from these 
ratios we also extract values for {f &/ fff) static ratios of couplings of excited states. The numerical 
results for the ratios are summarized in Table 2. 

6. Preliminary results for kinetic corrections 

Finally we present also some preliminary results for kinetic corrections obtained via the gen- 
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Figure 3: Ratios of meson decay constants (/bJ/b) static (left plot) and {fs'J Ib s ) static (right plot) as a func- 
tion of lattice spacing. The three quenched, continuum extrapolations correspond to (from left to right): a 
constant fit, a fit linear in a, a fit linear in a 2 . The green circle at a/ro « 0.3 corresponds to the adjacent 
quenched result, but with thin links for the static quark. 
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Figure 4: Preliminary results for the kinetic corrections. For the b quark mass we use m& w 4.2 GeV and 
assume Z « 1 for the renormalization constant. 

eralization of the variational method described in Section 2.3. Figure 4 shows first results for the 
G{\jm.b) corrections to the splitting 2S — IS, obtained from e' 11 ' and e^ 2,2 \ These are of course 
only rough estimates, where we have assumed that the renormalization factor of the operator QD 2 Q 
is close to unity, since we use a tadpole improved version of the operator D 2 . We also assume that 
the power-law divergent part of that operator is small (gt - ) 0) at this lattice spacing. As mass of 
the b quark we took a value of 4.2 GeV as an input. 



7. Summary 

In our studies, we use the variational method to extract not only the masses of ground and 
excited states but also the couplings. In addition, we generalize this method such that it allows also 
for extraction of three point functions of excited states. 

We successfully isolate excited static-light mesons via the variational method on a large num- 
ber of lattices, including 2S, 3S, IP, 2P, ID states. In general our results show quite good agreement 
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with the experimental values, where known. We perform also an extensive study for the static-light 
baryons. Also here we find a number of states. 

In addition to the spectrum, we try to determine decay constants of the ground and excited 
static-light mesons. So far we can only report ratios of the decay constants since the necessary 
renormalization constants are not determined yet for the particular actions which we use. 

As a next step we go beyond the static approximation by including kinetic corrections for the 
b quark. We treat these corrections as current insertions and develop a way to apply the varia- 
tional method to three -point functions involving excited states. First results have been presented, 
assuming a mild renormalization of the inserted current operator and using nib ~ 4.2 GeV. 
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